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ABSTRACT 

We present the bivariate brightness distribution (BBD) for the 2dF Galaxy Redshift 
Survey (2dFGRS) based on a preliminary subsample of 45,000 galaxies. The BBD 
is an extension of the galaxy luminosity function incorporating surface brightness 
information. It allows the measurement of the local luminosity density, js, and the 
galaxy luminosity and surface brightness distributions while accounting for surface 
brightness selection biases. 

The recovered 2dFGRS BBD shows a strong surface brightness-luminosity relation 
(Mb oc (2.4±J |)/x e ) providing a new constraint for galaxy formation models. In terms 
of the number-density we find that the peak of the galaxy population lies at Mb > 
— 16.0 mag. Within the well defined selection limits (—24 < Mb < —16.0 mag, 18.0 < 
/i e < 24.5 mag arcsec" 2 ) the contribution towards the luminosity-density is dominated 
by conventional giant galaxies (i.e. 90% of the luminosity-density is contained within 
—22.5 < M < —17.5, 18.0 < \i e < 23.0). The luminosity-density peak lies away from 
the selection boundaries implying that the 2dFGRS is complete in terms of sampling 
the local luminosity-density and that luminous low surface brightness galaxies are rare. 
The final value we derive for the local luminosity-density, inclusive of surface brightness 
corrections, is: jb — 2.49±0.20x lO 8 /iiooL Mpc~ 3 . Representative Schechter function 
parameters are: M* = -19.75 ± 0.05, <p* = 2.02 ± 0.02 x 10~ 2 and a = -1.09 ± 0.03. 
Finally we note that extending the conventional methodology to incorporate surface 
brightness selection effects has resulted in an increase in the luminosity-density of 
~ 37%. Hence surface brightness selection effects would appear to explain much of the 
discrepancy between previous estimates of the local luminosity density. 
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1 INTRODUCTION 



Of paramount importance in determining the mechanism(s) 
and epoch(s) of galaxy formation (as well as the local lumi- 
nosity density), is the accurate and detailed quantification 
of the local galaxy population. It represents the benchmark 
against which both environmental and evolutionary effects 
can be measured. Traditionally this research area originated 
with the all-sky photographic surveys coupled with a few 
handfuls of hard earned redshifts. Over the past decade this 
has been augmented by both CCD-based imaging surveys 
and multi-slit/fibre- fed spectroscopic surveys. From these 
data, a number of perplexing problems have arisen, most no- 
tably: the faint blue galaxy problem (Koo & Kron 1992; Ellis 
1997), the local normalisation problem (Maddox et al. 1990; 
Shanks 1990; Driver, Windhorst & Griffiths 1995; Marzke et 
al. 1998), the cosmological significance of low surface bright- 
ness galaxies (Disney 1976; McGaugh 1996; Sprayberry et 
al. 1996; Dalcanton et al. 1997; Impey & Bothun 1997) and 
dwarf galaxies (Babul & Rees 1992; Phillipps & Driver 1995; 
Loveday 1997; Babul & Ferguson 1996). These issues largely 
remain unresolved and arguably await an improved defini- 
tion of the local galaxy population (Driver 1999). 

Recent advancements in technology now allow for wide 
field-of-view CCD imaging surveys^ and bulk redshift sur- 
veys through purpose built multi-fibre spectrographs such as 
the common-user two-degree field (2dF) facility at the An- 
glo Australian Telescope (Taylor, Cannon & Parker 1998). 
The Sloan Digital Sky Survey elegantly combines these two 
facets (Margon 1999). 

The quantity and quality of data that is becoming avail- 
able allows not only the revision of earlier results but more 
fundamentally the opportunity to review and enhance the 
methodology with which the local galaxy population is rep- 
resented. For instance some criticism that might be levied at 
the current methodology — the representation of the space 
density of galaxies using the Schechter Luminosity function 
(Schechter 1976; Felten 1985; Binggeli, Sandage & Tam- 
mann 1988) — is that firstly, it assumes that galaxies are sin- 
gle parameter systems defined by their apparent magnitude 
alone, and secondly it describes the entire galaxy popula- 
tion by only three parameters; the characteristic luminosity 
L*, the normalisation of the characteristic luminosity <j>*, 
and the faint-end slope a. While it is desirable to represent 
the population with the minimum number of parameters, 
important information may lie in the nuances of detail. 

In particular, two recent areas of research suggest a 
greater diversity in the galaxy population than is allowed 
by the Schechter function form. Firstly, Marzke et al. (1994) 
and also Loveday (1997) report the indication of a change 
in the faint end slope at faint absolute magnitudes — a 
possible giant-dwarf transition — and this is also seen in a 
number of Abell clusters where it is easier to probe into the 
dwarf regime (e.g., Driver et al. 1994; De Propris et al. 1995; 
Driver, Couch & Phillipps 1998; Trentham 1998). Secondly a 
number of studies show that the three Schechter parameters, 
and in particular the faint-end slope, have a strong depen- 



* It is sobering to note that the largest published CCD-based 
imaging survey to date is 17.5 sq degrees to a central surface 
brightness of 25 V mag arcsec - 2 (Dalcanton et al. 1997) as com- 
pared to the all sky coverage of photographic media 
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Figure 1. Schechter luminosity functions from recent magnitude- 
limited redshift surveys (q.v. Table 1). The line becomes dotted 
outside the range of survey data. The range of values show the 
uncertainty in the LF which in turn filters through to the local 
measure of the mean luminosity density. 



dence upon: surface brightness limits (Sprayberry, Impey & 
Irwin 1996; Dalcanton 1998); colour (Lilly et al. 1996); spec- 
tral type (Folkes et al. 1999); optical morphology (Marzke 
et al. 1998), environment (Phillipps et al. 1998) and wave- 
length (Loveday 2000). It has been noted (Willmer 1997) 
that the choice of method for reconstructing the galaxy LF 
also contains some degree of bias. 

More fundamentally, evidence that the current method- 
ology might actually be flawed comes from comparing recent 
measurements of the galaxy luminosity function as shown in 
Fig. 1. The discrepancy between these surveys is significantly 
adrift from the quoted formal errors implying an unknown 
systematic error. The range of discrepancy can be quanti- 
fied as a factor of 2 at the L» (Af ~ —19.5) point rising 
to a factor of 10 at 0.01L, (M ~ -14.5). The impact of 
this variation is a factor of 3-4, for instance, in assessing the 
contribution of galaxies to the local baryon budget (e.g. Per- 
sic & Salucci 1992; Bristow & Phillipps 1994; and Fukugita, 
Hogan & Peebles 1998). 

This uncertainty is in addition to that introduced from 
the unanswered question of the space density of low surface 
brightness galaxies. The most recent attempt to quantify 
this is by O'Neil & Bothun (2000) — following on from Mc- 
Gaugh (1996), and in turn Disney (1976) — who conclude 
that the surface brightness function (SBF) of galaxies - 
the number density of galaxies in intervals of surface bright- 
ness — is of similar form to the luminosity function. Thus 
both the LF and SBF are described by a flat distribution 
with a cutoff at bright absolute magnitudes or high surface 
brightnesses. Taking the O'Neil result at face value, this im- 
plies a further error in measures of the local luminosity den- 
sity of 2-3 - i.e. the contribution to the luminosity- (and 
hence baryon-) density from galaxies is uncertain to a factor 
of ~ 10. However the significance of low surface brightness 
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galaxies depends upon their luminosity range and similarly 
the completeness of the LF relies on the surface brightness 
intervals over which each luminosity bin is valid. Both rep- 
resentations are incomplete unless the information is com- 
bined. This leads us to the conclusion that both the total 
flux and the manner in which this flux is distributed must 
be dealt with simultaneously. Several papers have been pub- 
lished which deal with either surface brightness distributions 
or Bivariate Brightness Distributions (Phillipps & Disney 
1986, Boyce & Phillipps 1994, Minchin 1999, Sodre & La- 
hav 1993, and Petrosian 1998). These are either theoretical, 
limited to cluster environments or have poor statistics due 
to the scarcity of good redshift data. 

Recently, Driver (1999) determined the first measure of 
the bivariate brightness distribution for field galaxies using 
Hubble Deep Field data (Williams et al. 1996) and capital- 
ising on photometric redshifts (Fernandez-Soto et al. 1998). 
The result, based on a volume limited sample of 47 galaxies, 
implied that giant low surface brightness galaxies were rare 
but that there exists a strong Luminosity- Surface Bright- 
ness relationship, similar to that seen in Virgo (Binggeli 
1993). The sense of the relationship implied that low surface 
brightness galaxies are preferentially of lower luminosity (i.e. 
dwarfs). If this is borne out it strongly tempers the conclu- 
sions of O'Neil & Bothun (2000). While the number of low 
surface brightness galaxies may be large, their luminosities 
are low, so their contribution to the local luminosity density, 
is also low, < 20% (Driver 1999). 

This paper attempts to bundle these complex issues 
onto a more intuitive platform by expanding the current rep- 
resentation of the local galaxy population to allow for: sur- 
face brightness detection effects, star-galaxy separation is- 
sues, surface brightness photometric corrections and cluster- 
ing effects. This is achieved by expanding the mono-variate 
luminosity function into a bivariate brightness distribution 
(BBD) where the additional dimension is surface brightness. 
The 2dFGRS allows us to do this for the first time by hav- 
ing a large enough database to separate galaxies in both 
magnitude and surface brightness without having too many 
problems with small number statistics. 

In §2 we discuss the revised methodology for measuring 
the space density of the local galaxy population, the local 
luminosity density and the contribution towards the baryon 
density in detail. In §3 we present the current 2dFGRS data 
(containing ~ 45, 000 galaxies or one fifth of the expected 
final tally). In §4 we correct for the light lost under the 
isophote and define our surface brightness measure. In §5 
we apply the methodology to construct the first statistically 
significant bivariate brightness distribution for field galaxies. 
The results for the number-density and luminosity- density 
are detailed in §6 and §7. In §8, we compare these results to 
other surveys. Finally we present our conclusions. 

Throughout we adopt H = 100kms _1 Mpc _1 and a 
standard fiat cosmology with zero cosmological constant 
(i.e. q — 0.5, A = 0). However we note that the results 
presented here are only weakly dependent on the cosmol- 
ogy- 



2 METHODOLOGY 

The luminosity density, j, is the total amount of flux emitted 
by all galaxies per Mpc 3 . When measured in the UV band 
it can be converted to a measure of the star-formation rate 
(see for example Lilly et al. 1996, Madau et al. 1998). When 
measured at longer wavelengths it can be combined with 
mass-to-light ratios to yield an approximate value for the 
contribution from galaxies towards the local matter density 
0,m — independent of H , only weakly cosmology dependent 
and not reliant on any specific theory of structure formation 
(see for example Carlberg et al. 1996; Fukugita, Hogan & 
Peebles 1998). The two main caveats are firstly the accuracy 
of jb (the luminosity density measured in the B-band) , and 
secondly the assumption of a ubiquitous mass-to-light ratio. 

2.1 Measuring j 

The luminosity density, j, is found by integrating the prod- 
uct of the number density $(L/L*) and the luminosity L 
with respect to luminosity. 

/•oo 

j= L $(L/L r )d(L/L r ) (1) 
Jo 

By convention, j is typically derived from a magnitude- 
limited redshift survey by determining the representa- 
tive Schechter parameters for a survey (e.g. Efstathiou et 
al. 1988) and then integrating the luminosity weighted 
Schechter function, where $(L/L») d(L/L,) is the Schechter 
function (Schechter 1976) given by: 

<b{L/L„)d{L/L,) = 0„ (L/L*) Q exp [-{L/L,)] d(L/L„) 

(2) 

and <f>*, I/*, and a are the three parameters which define the 
survey (referred to as the normalisation point, characteristic 
turn-over luminosity and faint-end slope parameter respec- 
tively). More simply if a survey is defined by these three 
parameters it follows that: 

j = <j>, L, r(a + 2) (3) 

Table 1 shows values for the luminosity density derived 
from a number of recent magnitude-limited redshift surveys 
(as indicated). The variation between the measurements of 
j from these surveys is ~ 2 and hence the uncertainty in the 
galaxy contribution to the mass budget is at best equally 
uncertain. This could be due to a number of factors, e.g., 
large scale-structure, selection biases, redshift errors, photo- 
metric errors or other incompleteness. In this paper we wish 
to explore the possibility of selection bias due to surface 
brightness considerations only. The principal motivation for 
this is that the LCRS (lower line on Fig. 1), which recov- 
ers the lowest j value, adopted a bright isophotal detection 
limit of )i r = 23 mag arcsec~ 2 , suggesting a dependence 
between the measured j and the surface brightness limit 
of the survey. Here we develop a method for calculating j 
which incorporates a number of corrections/considerations 
for surface brightness selection biases. In particular, a sur- 
face brightness dependent Malmquist correction, a surface 
brightness redshift completeness correction and an isophotal 
magnitude correction. We also correct for clustering. What 
is not included here, and will be pursued in a later paper, is 
the photometric accuracy, star-galaxy separation accuracy 
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and a detection correction specifically for the two degree 
field galaxy redshift survey (2dFGRS). 

Implementing these corrections requires re-formalising 
the path to j. Firstly, we replace the luminosity function 
representation of the local galaxy population by a bivariate 
brightness distribution (BBD). The bivariate brightness dis- 
tribution is the galaxy number density, $, as a function of 
absolute, total, B-band magnitude, Mb, and absolute, ef- 
fective surface brightness, fi e , i.e., <fr(M,n). To construct a 
BBD we need to convert the observed distribution to a num- 
ber density distribution taking into account the Malmquist 
bias and the redshift incompleteness correction, i.e., 



$(M, M ) 



V{M,n) 



W{M,v) 



(4) 



where: 



• 0(M, n) is the matrix of absolute magnitude, M, and 
absolute effective surface brightness, /i, for galaxies with red- 
shifts. 

• I(M, jj,) is the matrix of absolute magnitude, M, and 
absolute effective surface brightness, fi, for those galaxies 
for which redshifts were not obtained. 

• V(M, jj,) is the matrix which specifies the volume over 
which a galaxy with absolute magnitude, M, and abso- 
lute effective surface brightness, /i, can be seen (see also 
Phillipps, Davies & Disney 1990). 

• W(M, n) is the matrix that weights each bin to com- 
pensate for clustering. 

Deriving these matrices is discussed in detail later, j is then 
defined as: 



or in practice, 



M \± 

i = ^^L(M)<&(M,p) 

M M 



(5) 



^ = EE 10 ~ ' 4(M ~ Mo)$ ( M ^ ( 6 ) 

in units of L Mpc~ 3 where M 0B = 5.48. 

Our formalism has two key advantages over the tradi- 
tional luminosity function: Firstly, it adds the additional 
dimensionality of surface brightness allowing for surface 
brightness specific corrections. Secondly, it represents the 
galaxy population by a distribution rather than a function, 
thus requiring no fitting procedures or assumption of any 
underlying parametric form. 



arcsec -2 with a variation of ±0.11 mag arcsec -2 - Pimb- 
blet et al. 2000) . However APM b j magnitudes are found to 
vary from CCD bj magnitudes by 0.14±0.29mag (Metcalfe 
et al. 1995). Therefore the isophotal limit in APM bj magni- 
tudes is fium — 24.67 ± 0.30mag arcsec -2 . One pixel equals 
0.25 □". The minimum isophotal area found for galaxies in 
the APM catalogue is 35D". Star-galaxy separation was im- 
plemented as described in Maddox et al. (1990b). The final 
APM sample contains 3x 10 6 galaxies covering 15,000 square 
degrees - see Maddox et al. (1990a,b) for further details. 

The 2dFGRS input catalogue is a 2000 subregion 
of the APM catalogue (covering two continuous regions in 
the northern and southern Galactic caps plus random fields) 
with an extinction corrected magnitude limit of m = 19.45 
(Colless 1999). The input catalogue contains 250,000 galax- 
ies for which 81,895 have been observed using 2dF (as of 
November 1999). Each spectrum within the database has 
been examined by eye to check if the redshift is reliable. 
Redshifts are determined via cross-correlation with specified 
templates (see Folkes et al. 1999 for details). A brief test of 
the reliability of the 2dFGRS was achieved via a comparison 
between 1404 galaxies in common with the Las Campanas 
Redshift Survey (Lin et al. 1996) for which there were only 
8 mismatches, showing that 2dF redshifts are reliable. Of 
the 81,895 galaxies, 74,562 have a redshift, resulting in a 
redshift completeness of 91%. 

The survey comprises many overlapping two degree 
fields and many still have to be observed. Hence the abso- 
lute normalisation is tied to the full input galaxy catalogue 
which is known to contain 174.0 galaxies with m < 19.45 per 
square degree. Using a subsample of 44796 galaxies, cover- 
ing just the South Galactic Pole (SGP) region, this yields 
an effective coverage for this survey of 257 non-contiguous 
square degrees. 

Finally, for the purposes of this paper, we adopt a lower 
redshift limit of z = 0.015 to minimise the influence of pe- 
culiar velocities in the determination of absolute parameters 
and an upper redshift limit of z = 0.12. 

This upper limit of z = 0.12 was selected so as to max- 
imise the sample size yet minimise the error introduced by 
the isophotal corrections. At z = 0.12 the uncertainty in 
the isophotal correction (i^g), due to type uncertainty 
(see Appendix A), remains smaller than the photometric 
error (±0.2 mags). Note that the increase in the error in the 
isophotal correction is primarily because of the increase in 
the intrinsic isophotal limit due to a combination of surface 
brightness-dimming and the K-correction. 

The final sample is therefore pseudo volume- limited and 
contains 20,765 galaxies, with redshifts, selected from a par- 
ent sample of 45,000. 



3 THE DATA 

The data set presented here is based upon a sub-sample of 
the Automated Plate Measuring-machine galaxy catalogue 
(APM; Maddox et al. 1990a,b) for which spectra have been 
obtained using the two-degree field Facility (2dF) at the 
Anglo- Australian Telescope (AAT). 

The original APM catalogue contains bj magnitudes 
with random error Am = ±0.2 mag (Folkes et al. 1999) and 
isophotal areas Aiso- The isophotal area is defined as the 
number of pixels above a limiting isophote, /ii im — set at 
the 2cr level above the sky background (fiu m ~ 24.53 mag 



Note: all magnitude and surface brightnesses are in the APM 



bj filter. 



4 ISOPHOTAL CORRECTIONS 

The APM magnitudes have already been corrected assum- 
ing a Gaussian profile (see Maddox et al. 1990b for full de- 
tails) . This was aimed primarily at recovering the light lost 
due to the seeing and is crucial for compact objects. It is 
known to significantly underestimate the isophotal correc- 
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tion required for low surface brightness disks. Such systems 
typically exhibit exponential profiles with disks which can 
extend a substantial distance beyond the isophote, the most 
famous example being Malin 1 (Bothun et al. 1987). Once 
thought of as a Virgo dwarf this system remains the most 
luminous field galaxy known. 

To compliment the Gaussian correction (required for 
compact objects but ineffectual for extended sources) we 
introduce an additional correction (ineffectual for compact 
sources but suitable for extended disks). This correction as- 
sumes all objects can be represented by a pure exponential 
surface brightness profile extending from the core outwards. 
In this case the surface brightness profile is simply: 

E(r) = E cxp(-r/a) (7) 

or, 

fi(r) = fi ± 1.086(r/a) (8) 

Where E D is the central surface brightness in Wm~ 2 
arcsec -2 , a is the scale length of the galaxy in arcsecs and 
r the radius in arcsecs. /i is the central surface brightness 
in mag arcsec -2 . 

Under this assumption a galaxy's observed isophotal 
luminosity is the integrated radial profile out to ri SO , 

f r iso 

Uso — 2tv / E D exp(— r/a)rdr (9) 
Jo 

which can be expressed in magnitudes as: 

m lao = ^" pp - 2.51og 10 {2 7r[a 2 - ct(a + r lso ) exp(-r lso /a)]} 

(10) 

(here fj,„ pp denotes the apparent surface brightness un- 
corrected for redshift.) \iu-m, the detection/photometry 
isophote, can be expressed as: 

VUm = n7 P + 1.086(r iso /a) (11) 

As rriiso, ri 30 and uu m are directly measurable quantities, 
equations ( [lo] & jll] ) can be solved numerically. The total 
magnitude is then given by: 

|>OC 

l tot = 2TT T, exp(-r/a)rdr = 27rE D a 2 (12) 
Jo 

or, 

m tot =^-2.51og 10 (2™ 2 ) (13) 

From this description an extrapolated central surface 
brightness can be deduced numerically from the specified 
isophotal area and isophotal magnitude (after the seeing 
correction). Note that this prescription ignores the possible 
presence of a bulge, opacity, and inclination leading to an un- 
derestimate of the isophotal correction. This is unavoidable 
as the data quality is insufficient to establish bulge-to-disk 
ratios. To verify the impact of this we explore the accuracy 
of the isophotal correction for a variety of galaxy types in 
Appendix A. The tests show that the isophotal correction 
is a significant improvement over the isophotal magnitudes 
for all types - apart for ellipticals where the introduced er- 
ror is negligible compared to the photometric error - and 
a dramatic improvement for low surface brightness systems. 
The final magnitudes, after isophotal correction, now lie well 
within the quoted error of ±0.2 mags for both high- and low- 
surface brightness galaxies. 



4.1 The Effective Surface Brightness 

Most results cited in the literature use the central surface 
brightness or the effective surface brightness. The central 
surface brightness, as described above, is the extrapolated 
surface brightness at the core under the assumption of a per- 
fect exponential disk. The effective surface brightness is the 
mean surface brightness within the half-light radius. The 
conversion between the measures is relatively straightfor- 
ward and described as follows: 

li = 7r E a 2 = 2ir E a 2 [l — (1 + r e /a) exp(— r e /a)] 

(14) 

which can be solved numerically to get 

r e = 1.678a (15) 
The effective surface brightness is now given by: 

li a / p = ii7 v + 2.51og 10 [(r e /a) 2 ] = tf* + 1.124 

(16) 

Hence from the isophotal magnitudes and areas we can de- 
rive the total magnitude and effective surface brightness 
(quantities which are now independent of the isophotal de- 
tection threshold). We chose to work with effective surface 
brightness as it can, at some later stage, be measured di- 
rectly from higher quality CCD data. Note that these surface 
brightness measures are all apparent rather than intrinsic, 
however this is not important as although surface brightness 
is distance dependent the isophotal correction is not (this is 
because both fj,u m and /i e vary with redshift in the same 
way). 

Fig. 2 shows the final 2dFGRS sample (i.e. after isopho- 
tal correction) for those galaxies with (upper panel) and 
without (lower panel) redshifts. The galaxies are plotted ac- 
cording to their apparent total magnitude and apparent ef- 
fective surface brightness. The curved boundary at the faint 
end of both plots is due to the isophotal corrections which 
are strongly dependent on fj, e for a constant m. As /x e — 1.124 
tends towards (J,u m , the isophotal limit, the correction tends 
towards infinity, making it impossible to see galaxies close 
to fiu m . The average isophotal correction is 0.33 mag (for 
Hum = 24.67 mag arcsec -2 ). 

The observed mean magnitude and observed mean ef- 
fective surface brightness for those galaxies with and without 
redshifts are: 18.06 ±0.01 mag & 22.66 ±0.01 mag arcsec -2 
and 18.54 ± 0.01 mag & 23.17 ± 0.01 mag arcsec -2 , respec- 
tively. These figures imply that galaxies closer to the detec- 
tion limits are preferentially under-sampled. 



5 CONSTRUCTING THE BBD 

We now apply the methodology described in section 2 to 
derive the BBD from our data set. This requires constructing 
the four matrices, 0(M,/x), I(M,fj), V(M, fj,) and W(M,n). 

5.1 Deriving 0{M, fj,) 

For those galaxies with redshifts, we can obtain their ab- 
solute magnitude and absolute effective surface bright- 
ness assuming a cosmological framework and a global K- 
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Figure 2. Galaxies with (upper) and without (lower) redshifts 
for the current 2dFGRS sample plotted according to their appar- 
ent extinction corrected total magnitude and apparent effective 
surface brightness. 



Figure 3. Galaxies from the 2dFGRS with redshifts, plotted 
in absolute magnitude and absolute effective surface brightness 
space. The shaded region denotes the regions where less than 10 4 
Mpc 3 are surveyed and is based on visibility theory as described 
in Phillipps, Davies & Disney (1990). The 3 curves represent the 
volumes of 10 4 , 10 5 and 5 X 10 5 Mpc 3 . 



correction^ K(z)=2.5z (Driver et al. 1994). The conversions 
from observed to absolute parameters are given by: 

2c 



M 



and, 



5 log! 



■(1 +«)(!- (l + ^r°' b )]-25 



Me = Me 



10Iog 10 (l + *)-*(«) 



K(z) 
(17) 



(18) 



Here H is the Hubble constant, c is the speed of light, 
Me PP is the apparent effective surface brightness and /i e is 
the absolute effective surface brightness. The M, derived by 
Eqn is a total absolute magnitude since the correction 
has been made for the light below the isophote. 

Fig. 3 shows the upper panel of Fig. 2 with the axes 
converted to absolute parameter space using the conver- 
sions shown above. Naturally galaxies in different regions 
are seen over differing volumes, because of Malmquist bias, 
hence it is not yet valid to compare the relative numbers. 
However it is possible to define lines of constant volume 
as shown on Fig. 3 (dotted lines). These lines are derived 
from visibility theory (Phillipps, Davies & Disney 1990) 
and they delineate the region of the BBD plane where 
galaxies can be seen over various volumes. The shaded re- 
gion shows the region where the volume is less than 10 4 



t Individual K-corrections will be derived from the data; however 
this has not yet been implemented. 
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Mpc 3 and hence where we are insensitive to galaxy densi- 
ties of < 10 -2 galaxies/Mpc -3 mag -1 (mag arcsec -2 ) -1 . The 
equations used to calculate the lines are laid out in Ap- 
pendix B. We show a V = 5 x 10 5 Mpc 3 line rather than 
V = 10 6 Mpc 3 because the z = 0.12 limit is at a volume less 
than V = 10 6 Mpc 3 . The parameters used in the visibility 
calculations are: fxiim = 24.67 mag arcsec -2 ; 9 min = 7.2", 
m ax = 200.0", mbright = 14.00 mag; m faint = 19.45 mag; 
z-min = 0.015; and z max = 0.12. The clear space between the 
data and the selection boundary at bright absolute magni- 
tudes implies that although the 2dFGRS data set samples 
a sufficiently representative volume (V > 10,000 Mpc 3 ), 
galaxies only exist over a restricted region of this observed 
BBD. Fainter than M — —16.5 mags the volume is insuf- 
ficient to sample populations with a space density of 10 -2 
Mpc -3 mag -1 (mag arcsec -2 )" 1 or less. 

Fig. 4 shows the data of Fig. 3 binned into /i e and M in- ^ 
tervals to produce the matrix 0(M, fx) (see EqnQ). The bins 
are 0.5 magx0.5 mag arcsec -2 and start from — 24.0mag and 
a central surface brightness of 19.0mag arcsec -2 , effective 
surface brightness of 20.12mag arcsec -2 . The total number 
of bins is 200 in a uniform 20 x 10 array. Fig. 4 represents 
the observed distribution of galaxies and shows a strong peak 
close to the typical M* value seen in earlier surveys (see ref- 
erences listed in Table 1). 

5.2 Deriving I(M, y) 

Not all galaxies targeted by the 2dFGRS have a measured 
redshift. This may be due to lack of spectral features, se- 
lection biases or a misplaced/defunct fibre. One method to 
correct for these missing galaxies is to assume that they have 
the same observed BBD as those galaxies for which redshifts 
have been obtained. One can then simply scale up all bins 
by this known incompleteness. 

However, from Fig. 2 and §1 & 3 we noted that the 
incompleteness is a function of both the apparent magni- 
tude and the apparent surface brightness. There is no reli- 
able way of converting these values to absolute values with- 
out redshifts and to obtain an incompleteness correction, 
I(M,/j,), some assumption must be made. Here we assume 
that a galaxy of unknown redshift with apparent magni- 
tude, m, and apparent effective surface brightness, has 
a range of possible BBD bins that can be statistically rep- 
resented by the BBD distribution of galaxies with m ± Am 
and fj, e ± A fie- The underlying assumption is that galaxies 
with and without redshifts with similar observed m and fi 
have similar redshift distributions. I.e., the detectability of 
a galaxy is primarily dependent on its apparent magnitude 
and apparent surface brightness. [While these factors are ob- 
viously crucial one could also argue that additional factors, 
not incorporated here, such as the predominance of spec- 
tral features are also important such that the true probabil- 
ity distribution for the missing galaxies could be somewhat 
skewed from that derived in Fig. 5.] 

Hence for each galaxy, without a redshift, we select 
those galaxies, with redshifts, within O.lmag and O.lmag 
arcsec -2 and determine their collective observed BBD dis- 
tribution. This is achieved using all 44796 galaxies in the 
SGP sample, as galaxies, without redshifts, are not limited 
to z < 0.12. This distribution is then normalised to unity to 
generate a probability distribution for this galaxy. This is re- 
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Figure 4. The observed distribution of the 2dFGRS data set 
mapped onto the BBD, prior to volume and incompleteness cor- 
rection. The contour lines are set at 25, 100, 250, 500, 750, 1000, 
1250, 1500 and 1750 galaxies per bin. The minimum number of 
galaxies in a bin is 25. 

peated for every galaxy without a redshift. Fig. 5, shows the 
probability distribution for a galaxy with m = 17.839 and 
He = 22.314. There are 317 galaxies with redshifts within 
O.lmag and O.lmag arcsec -2 and combined they have a dis- 
tribution ranging from M = —15.44 to M = —21.76 and 
H e = 22.26 to He = 21.06. 

To generate the matrix, I(M, p), the probability dis- 
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Figure 5. The BBD distribution of galaxies with m = 17.84±0.1 
and iil vv = 22.31 ± 0.1. The contour lines are set at 0.01, 0.05, 
0.10, 0.15 and 0.20 chance of the galaxy being in that bin. 



tributions for each of the galaxies without redshifts, (nor- 
malised to unity) are summed together to give the total 
distribution of those galaxies without redshifts. However 
these galaxies could have the full range of redshifts that 
each galaxy can be detected over, not a range limited to 
z — 0.12. Therefore the number distribution is weighted by 
the fraction of galaxies within each bin with z < 0.12. 

Fig. 6 shows I(M,fx) which should be compared to Fig. 
4 (0(M, //)). The distribution of Fig 6 appears broader in- 



dicating that the missing galaxies are not random but that 
they are predominantly low luminosity, low surface bright- 
ness systems. This is illustrated in Fig. 7, which shows the 
ratio of I(M, fi) to 0(M, jj.) for the bins containing more 
than 25 galaxies with redshifts. From Fig. 7 we see that the 
trend is for the ratio to increase towards the low surface 
brightness regime. There is no significant trend in absolute 
magnitude. Finally we note that the although the incom- 
pleteness correction does increase the population in some 
bins by as much as 50% we shall see in §6 & 7 that the con- 
tribution from these additional systems towards the overall 
luminosity density is negligible. 

5.3 Deriving V(M, y) 

To convert the number of observed galaxies to number 
density per Mpc 3 a Malmquist correction is required, i.e., 
V(M,y). This matrix reflects the volumes over which each 
M, fj, bin can be observed. One option is to use visibility the- 
ory as prescribed by Phillipps, Davies & Disney (1990) - and 
used to construct the constant volume line on Fig. 3. While 
visibility is clearly a step in the right direction, and prefer- 
able to applying a magnitude-only dependent correction, its 
limitation is that it assumes idealised galaxy profiles (i.e. 
it neglects the bulge component, seeing, star-galaxy sepa- 
ration and other complications). Ideally one would like to 
extract the volume information from the data itself and this 
is possible by using a 1/VWx type prescription, i.e., within 
each 0(M, y) bin, the maximum redshift at which a galaxy 
is seen is determined and the volume derived from this red- 
shift. The advantages of using the data set rather than the- 
ory is that it naturally incorporates all redshift dependent 
selection biases. However the maximum redshift is suscep- 
tible to scattering from higher-visibility bins. An improved 
version is therefore to use the 90*' 1 percentile redshift and 
to adjust 0(M,n) and I(M,y) accordingly. 

Although this requires rejecting 10% of the data it has 
two distinct advantages. Firstly it ensures that the redshift 
distribution in each bin has a sharp cutoff (as opposed to 
a distribution which peters out). Secondly it uses the entire 
dataset as opposed to the maximum redshift only. In the 
case where the 90 th percentile is not exact we take the galaxy 
nearest. Using these redshifts the volume can be calculated 
independently for each bin assuming an Einstein-de Sitter 
cosmology as follows: 



V{M,n) 



T/ Z 90( M >M) 
V z min (M, f i) 



where, 



v. 



z 90 (M,ti) 
,(M, M ) 



z 90 (M, M ) 



»(M, M ) 



H (1 + Z )»' 5 



dz 



(19) 



(20) 



and a is the solid angle in steradians on the sky. di is the 
luminosity distance to the galaxy. z m in is the minimum dis- 
tance over which a galaxy can be detected. z m i„ is calculated 
from visibility theory (Phillipps, Davies & Disney 1990, see 
also Appendix B) adopting values for the maximum magni- 
tude and maximum size of tub = 14.00 mag and 6 — 200" 
respectively. 

Fig. 8 shows the matrix ^ . Note that squares con- 
taining fewer than 25 galaxies are not shaded. This matrix 
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Figure 6. The bivariate number distribution of those galaxies 
without redshifts, i.e., I(M,fi). The contour lines are set at 10, 
25, 50, 75, 100, 125 and 150 galaxies bin" 1 . 

is flat-bottomed due to the cutoff at z = 0.12. Fig. 8 shows a 
strong dependency upon magnitude (i.e. classical Malmquist 
bias as expected) and also upon surface brightness. This 
surface brightness dependency is particularly strong near 
the 10 4 Mpc 3 volume limit, as the data becomes sparse. In- 
side this volume limit the contour lines generally mimic the 
curve of the visibility-derived volume boundary. This sug- 
gests that visibility theory provides a good description of 
the combined volume dependency. The sharp cutoff along 
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Figure 7. The ratio of galaxies without redshifts to galaxies with 
redshifts. The contour lines arc set at 0.01, 0.05, 0.10, 0.15, 0.20, 
0.25, 0.30, 0.35, 0.40, 0.45 and 0.50. 

the high surface brightness edge may be real but could also 
be a manifestation of the complex star-galaxy separation al- 
gorithm (see Maddox 1990a). Given that a galaxy seen over a 
larger distance appears more compact and that local dwarfs 
have smaller scale lengths than giants (cf. Mateo 1998), this 
seems reasonable. We will investigate this further through 
high-resolution imaging. 

The main point to take away from Fig. 8 is that the visi- 
bility surface of the 2dFGRS input catalogue is complex and 
dependent on both M and jj, e (although predominantly M). 
Any methodology which ignores surface brightness informa- 
tion and implements a volume-bias correction in luminos- 
ity only, is implicitly assuming uniform visibility in surface 
brightness. The 2dFGRS data clearly show that this is not 
the case. 

5.4 Clustering 

Structure is seen on the largest measurable scales (e.g. de 
Lapparent, Cellar & Huchra 1986). To determine whether 
the effects of clustering are significant we constructed a 
radial density profile, as shown in Fig. 9. This was de- 
rived from those bins for which more than 100 galaxies 
are seen over the whole range from 0.015 < z < 0.12 (i.e. 
-19.75 < M < -18.75 and 21.1 < ^ e < 22.6). Those galax- 
ies which are brighter cannot be seen at z < 0.015 due to the 
bright magnitude cut a,t m — 14.00 and would therefore bias 
the number density towards the bright end. For these high- 
visibility galaxies we calculated their number-density {<f>) in 
equal volume intervals of 5.0 x 10 3 Mpc 3 , from z = 0.0185 
to z = 0.12. Fig. 9 shows that clustering is severe with what 
appears to be a large local void around z = 0.04 and walls 
at z = 0.06 and z = 0.11. The ESO Slice Project (ESP) 
survey (Zucca et al. 1997) whose line-of-sight (RA ~ OOh, 
5 ~ —40) is just outside the 2dF SGP region, measures an 
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Figure 9. The clustering map. This shows the number density of 
giant galaxies (-19.75 < M < -18.75 and 21.1 < fJ. e < 22.6 as a 
function of redshift. The points are spaced equally in volume at 
intervals of 5000 h 3 Mpc 3 starting at 10000 h 3 Mpc 3 . 



visibility galaxies by taking the ratio of the number-density 
of high- visibility galaxies over the full redshift range divided 
by the number-density of high-visibility galaxies over the 
redshift range of each bin, i.e.: 



W(M,fi) = 



4>{High - 



• \z=0.12 

- «*s) 2=0 . 015 



<j)(High — vis) 



zmax(M \[t) 
zmin(M 



(21) 



Absolute Magnitude / mag 



This weighting matrix is shown in Fig. 10. The impli- 
cation is that the number-density of low-luminosity systems 
will be amplified by almost a factor of 1.5, to correct for 
the apparent presence of a large local void along the SGP 
region - indicated by the vertical ridge at M = —17 on Fig. 
10. Once again this implicitly assumes that the clustering of 
dwarf and giant systems is correlated. 



Figure 8. Plot of 1/volume in each bin. The contour lines are 
set at 1.0 x 10~ 7 , 1.0 x 10~ 6 , 2.0 x 10~ 6 , 1.0 x 10~ 5 , 2.0 x 10" 5 , 
4.0 x 10" 5 , 6.0 x 10~ 5 , 8.0 x 10" 5 and 1.0 x 10" 4 h" 3 Mpc" 3 . 

under-density at < 140/i _1 Mpc [z « 0.045) and an over- 
density at z ~ 0.1. The structure that they see closely re- 
sembles the structure that we see. A reliable measure of the 
BBD needs to correct for this clustering-bias. Here we adopt 
a strategy which implicitly assumes, firstly, that clustering 
is independent of either M or y, and secondly, that evolu- 
tionary processes to z = 0.12 are negligible. 

On the basis of these caveats we constructed a weight- 
ing matrix, W{M,y). This was determined from the high- 



6 THE 2DFGRS BBD 

Finally we can combine the four matrices, 0(M, y), I(M, y), 
V(M, fi) and W(M, y) (see Eqn |), to generate the 2dFGRS 
bivariate brightness distribution, as shown in Fig. 11. This 
depicts the underlying local galaxy number-density distri- 
bution inclusive of surface brightness selection effects. Only 
those bins which are based upon 25 or more galaxies are 
shown. Note that by summing the BBD along the surface 
brightness axis one recovers the luminosity distribution of 
galaxies. By summing along the magnitude axis one obtains 
the surface brightness distribution of galaxies (see §8). 

Fig. 12 shows the errors in the BBD. These were initially 
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Figure 10. The Weighting Map. The contours are at 0.8, 0.9, 
1.0, 1.1, 1.2, 1.3, 1.4 and 1.5. 



determined via Monte-Carlo simulations, assuming a Gaus- 
sian error distribution of ±0.2 mag, in the APM mag nitudes . 
This showed that the errors were proportional to \J (1/N), 

and since (1/N) is much faster to calculate, this is the 
result that is used throughout the calculations. These errors 
were then combined in quadrature with the additional er- 
ror in the volume estimate, assuming Poisson statistics. The 
total error is given by: 
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Figure 11. The 2dFGRS Bivariate Brightness Distribution. The 
contour lines are set at 1.0 x 10~ 7 , 1.0 x 10~ 3 , 2.5 x 10~ 3 , 5.0 x 
10- 3 , 7.5 x 10- 3 , 1.0 x 10~ 2 , 1.25 x lO" 2 , 1.5 x lO" 2 , 1.75 x 10~ 2 , 
2.0 X 10~ 2 , and 2.25 X 10~ 2 galaxies Mpc~ 3 bin" 1 . The thick 
lines represent the selection boundaries calculated from visibility 
theory. 



a * = \l ( N tot (M^) )+( wk^ ] (22) 

The errors become significant (> 20%), when M > —16 
and around the boundaries of the BBD shown in Fig. 11. 
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Figure 12. The errors in the BBD. The contour lines are set 
at 1.0 x 10~ 7 , 1.0 x 10~ 4 , 5.0 x 10~ 4 , 1.0 x 10~ 3 , 1.5 x 10~ 3 , 
2.0 x 10~ 3 , 2.5 x 10~ 3 , 3.0 x 10~ 3 , 3.5 x 10~ 3 and 4.0 x 10~ 3 
galaxies Mpc~ 3 bin -1 . 



The data and associated errors are tabulated in Table 2. 
From Figs 11 & 12 we note the following: 



6.1 A luminosity-surface brightness relation 

The BBD shows evidence of a luminosity-surface brightness 
relation similar to that seen in Virgo, (Mb oc 1.6/ito, Binggeli 
1993), in the Hubble Deep Field, (Mfasow oc 1.5/i e , Driver 
1999) and in Sdm galaxies (M B oc 2.02 ± 0.16/U e de Jong 
& Lacey 2000). A formal fit to the 2dFGRS data yields 
Mb oc (2.4±J;|)/i e . While confirming the general trend this 
result appears significantly steeper than the Virgo and HDF 
results. Both the Virgo and HDF results are based on lower 
luminosity systems hence this might be indicative of a sec- 
ond order dependency of the relation upon luminosity. Alter- 
natively it may reflect slight differences in the data/analysis 
as neither the Virgo nor HDF data include isophotal correc- 
tions whereas the 2dFGRS data is more susceptible to at- 
mospheric seeing. The gradient is slightly steeper than the 
de Jong & Lacey result, but is well within the errors. 

The presence of a luminosity-surface brightness relation 
highlights concerns over the completeness of galaxy surveys, 
as surveys with bright isophotal limits will preferentially ex- 
clude dwarf systems leading to an underestimate of their 
space-densities and variations such as those seen in Fig. 1. 

The confirmation of this luminosity-surface brightness 
relation within such an extensive dataset is an important 
step forward and any credible model of galaxy formation 
must now be required to reproduce this relation. 



6.2 A Dearth of Luminous Low surface brightness 
galaxies 

Within each magnitude interval there appears to be a pre- 
ferred range in surface brightness over which galaxies may 
exist. While the high-surface brightness limit may be due, 
in part or whole, to star-galaxy separation and/or fibre- 
positioning accuracy the low surface brightness limit ap- 
pears real. This cannot be a selection limit as one requires 
a mechanism which hides luminous LSBGs yet allows dwarf 
galaxies of similar surface brightness to be detected within 
the same volume. The implications are that these galaxy 
types (luminous-LSBGs) are rare with densities less than 
10~ 4 galaxies Mpc -3 . This result is important as it directly 
addresses the issues raised in the introduction and implies 
that existing surveys have not missed large populations of 
luminous low surface brightness galaxies. Perhaps more im- 
portantly it confirms that the 2dFGRS is complete for giant 
galaxies and that the postulate that the Universe might be 
dominated by luminous LSBGs (Disney 1976) is ruled out. 

One caveat however is that luminous-LSBGs could be 
masquerading as dwarfs. For example consider the case of 
Malin 1 (Bothun et al. 1987) which has a huge extended 
disk (55 kpc scale-length) of very low surface brightness 
(/i = 26.5). This system is actually readily detectable be- 
cause of its high surface brightness active core, however 
within the 2dFGRS limits it would have been miss-classified 
as a dwarf system with M = —17.9, /i e = 21.8. Hence Figs 
11 & 12 rule out luminous disk systems only. To determine 
whether objects such as Malin 1 are hidden amongst the 
dwarf population will require either ultra-deep CCD imag- 
ing or cross-correlation with Hi-surveys which would exhibit 
very high HI mass-to-light ratios for such systems. 



6.3 The rising dwarf population 

The galaxy population shows a steady increase in number- 
density with decreasing luminosity. This continues to the 
survey limits at M = —16, whereupon the volume limit and 
surface brightness selection effects impinge upon our sample. 
The expectation is that the distribution continues to rise 
and hence the location of the peak in the number-density 
distribution remains unknown. However we do note that the 
increase seen within our selection limits is insufficient for 
the dwarf population to dominate the luminosity-density as 
shown in the next section. 

Perhaps more surprising is the lack of sub-structure in- 
dicating either a continuity between the giant and dwarf 
populations or that any sub-structure is erased by the ran- 
dom errors. The former case is strongly indicative of a hier- 
archical merger scenario for galaxy formation which one ex- 
pects to lead towards a smooth number-density distribution 
between the dwarf and giant systems (White & Rees, 1978). 
This is contrary to the change in the luminosity distribution 
of galaxies seen in cluster environments (e.g. Smith, Driver, 
& Phillipps 1997). In a later paper we intend to explore the 
dependency of the BBD upon environment. 
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7 THE VALUE OF J B AND Sl A 



The luminosity density matrix, j(M,fj.) is constructed from 
L <E>(A/, fi) in units of LqMpc -3 and is shown as Fig. 13. 
The distribution is strongly peaked close to the conventional 
M* parameter derived in previous surveys (see Table 1). 
The peak lies at M = —19.5 mag and /j, e = 22.12 mag 
arcsec - . The final value obtained is js = (2.49 ± 0.20) x 
10 8 /iioo£©Mpc~ 3 . The sharp peak sits firmly in the centre of 
our observable region of the 2dFGRS BBD and drops rapidly 
off on all sides towards the 2dFGRS BBD boundaries. This 
implies that while the 2dFGRS does not survey the entire 
parameter space of the known BBD it does effectively con- 
tain the full galaxy contribution to the local luminosity den- 
sity. Redoing the calculations using galaxies with redshifts 
also results in j B = (2.49 ± 0.20) x IO^iooI/qMpc -3 . This 
demonstrates that there is no dependency of the results upon 
the assumption made for the distribution of galaxies without v 
redshifts. 

We note that j derived via a direct — ^ — estimate, with- 
out any surface brightness or clustering corrections, (Eqn^) 
gives a value of j = 1.82 ± 0.07 x lO 8 /iiooZ/0Mpc" 3 . Includ- 
ing the isophotal magnitude correction only leads to a value 
of j = 2.28 ± 0.09 x 10 8 /iioo£©Mpc" 3 . Hence a more de- 
tailed analysis leads to a 36.8% increase in j, of this 36.8%, 
25.3% is due to the isophotal correction, 10.4% is due to the 
Malmquist bias correction and 1.1% is due to the clustering 
correction. 

The final value agrees well with that obtained from the 
recent ESO Slice Project (Zucca et al. 1997). The method 
that they used corrects for clustering, but not surface bright- 
ness — although their photometry is based on aperture 
rather than isophotal magnitudes. However it is worth point- 
ing out that j is sensitive to the exact value of fj,u m . The 
quoted error in fj,u m is ±0.3 (for plate-to-plate variations, 
Metcalfe et al 1995; Pimblett et al 2000). Table 3 shows a 
summary of results when repeating the entire analysis using 
the upper and lower error limits. It therefore seems likely 
that a combination of surface brightness biases and large 
scale structure can indeed lead to the type of variations seen 
in Table 1. 

Finally following the method of Carlberg, Yee & Elling- 
son (1997) we can obtain a crude ball-park figure for the to- 
tal local mass density by adopting a universal mass-to-light 
ratio based on that observed in clusters. While this method 
neglects biasing (White, Tully & Davis, 1988) it does pro- 
vide a useful crude upper limit to the mass density. From 



Carlberg et al. (1997) we find: 



= 289 ± 50hi 



Assuming a mean colour of (B — R) = 



index of 1.17 this converts to: 



1.1 and a solar colour 

l b ' = 271±47/uoo^. 
Multiplying the luminosity density by the mass-to- light ratio 
yields a value for the local mass-density of Qm ~ 0.24. We 
note that this is consistent with the current constraints from 
the combination of Sn la results with the recent Boomerang 
and Maxima-1 results (Balbi et al. 2000; de Bernardis et al. 
2000). 



8 COMPARISONS WITH OTHER SURVEYS 

As this work represents the first detailed measure of the 
field BBD there is no previous work with which to com- 




-20 -18 

Absolute Magnitude / mag 

Figure 13. Plot of the Luminosity Density Distribution. The 
contour lines are set at 100, l.OxlO 6 , 5.0xl0 6 , l.OxlO 7 , 2.0xl0 7 , 
3.0 x 10 7 , 4.0 x 10 7 , 5.0 x 10 7 , and 6.0 x 10 7 L Q Mpc" 3 bin" 1 . 



pare. However as mentioned earlier it is trivial to convert the 
BBD into either a luminosity distribution and/or a surface 
brightness distribution both of which have been determined 
by numerous groups. This is achieved by summing across 
either luminosity or surface brightness intervals. For those 
bins containing fewer than 25 galaxies for which a volume- 
and clustering- bias correction were not obtained we use 
the volume-bias correction from the nearest bin with 25 or 
more galaxies. This will lead to a slight underestimate in 
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the number-densities, however as the number-density peak 
is well defined this effect is negligible. 



8.1 The luminosity distribution 

Fig. 14 shows a compendium of luminosity function mea- 
sures (see Table 1 and Fig. 1). Superimposed on the previ- 
ous Schechter function fits (dotted lines) are the results from 
the 2dFGRS. Three results are shown, the luminosity dis- 
tribution neglecting surface brightness and clustering (short 
dashed line), the luminosity distribution inclusive of sur- 
face brightness corrections (long dashed line) and finally the 
luminosity distribution inclusive of surface brightness and 
clustering corrections (solid line) . The inset shows the 2a x 2 
error eclipse for the latter case, yielding Schechter function 
parameters of M* = -19.75 ± 0.05, a = -1.09 ± 0.03 and 
<t>* = 0.0202 ± 0.0002 in close agreement with the ESO Slice 
Project (ESP, see Table 1). 

Comparing the 2dFGRS result with other surveys sug- 
gests that the main effect of the surface brightness correction 
is to shift M* brightwards by the mean isophotal correction 
of 0.33 mag and to increase the number density by a factor 
of 1.2. The clustering correction has little effect at bright 
magnitudes, almost by definition, but significantly amplifies 
the dwarf population. This is a direct consequence of an ap- 
parent local void at z=0.04 along the full range of the 2dF 
SGP region. 

We note that the combination of these two corrections 
mimics closely the discrepancies between various Schechter 
function fits. For example the shallower SSRS2, APM and 
Durham/UKST surveys are biased at all magnitudes by the 
apparent local void while the deeper ESP and Autofib sur- 
veys (which employ clustering independent methods) probe 
beyond the void. Similarly because of the luminosity- surface 
brightness relation those surveys which probe to lower sur- 
face brightnesses will find a higher dwarf-to-giant ratio. 

The shaded region shows the point at which our data 
starts to become highly uncertain because of the small vol- 
ume surveyed (V < 10, 000 Mpc 3 ). As this is the first statis- 
tically significant investigation into the bivariate brightness 
distribution we can conclude that as yet no survey contains 
any direct census of the space density of Mb > —16.0 mag 
galaxies. 
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Figure 14. The solid line shows the final luminosity function, 
with all corrections taken into account. The long dashed line 
shows the luminosity function with surface brightness corrections 
only, and the short dashed line shows the luminosity function ig- 
noring surface brightness and clustering corrections. Also shown, 
as dotted lines, are the luminosity functions from Fig. 1. The 
shaded region denotes the limit of reliability based on visibility 
theory. 



the range 22 < [i < 25 albeit with substantial scatter. Fig. 
15 reproduces Fig. 1. of O'Neil & Bothun but now includes 
the final 2dFGRS results. In order to compare our results 
directly with the O'Neil & Bothun data we assumed a mean 
bulge-to-total ratio of 0.4 (see Kent 1985) resulting in a uni- 
form offset of 0.55 mag arcsec -2 . 

The 2dFGRS data is substantially broader than the 
Freeman distribution and appears to agree well with the 
compendium of data summarised in O'Neil & Bothun 
(2000). From visibility theory (see Fig. 3) we are com- 
plete (i.e. the volume observed is greater than 10 4 Mpc 3 
and therefore statistically representative) in central surface 
brightness from 18.0 < pi < 23 for M < —16. Assuming 
that the luminosity-surface brightness relation continues as 
reported in Driver (1999) and Driver & Cross (2000) the 
expectation is that the surface brightness distribution will 
steepen as galaxies with lower luminosities are included. 



8.2 The surface brightness distribution 

There have been a few Surface Brightness Functions pub- 
lished over the years. The first was the Freeman (1970) re- 
sult, which showed a Gaussian distribution with fi — 21.65 
and o-jj = 0.3. Since then, however, many galaxies have been 
found at greater than 10cr from the mean. The probability of 
a galaxy occurring at lOcr or greater is w 1 x 10 -20 . The total 
number of galaxies in the universe is in the range of 10 11 — 
10 12 , so the Freeman SBF must be underestimating the LS- 
BGs. The distribution seen by Freeman is almost certainly 
due to the relatively bright isophotal detection threshold 
that the observations were taken at, around 22 - 23 mag 
arcsec -2 . A more recent measure of the SBF comes from 
O'Neil & Bothun (2000) who also show a compendium of 
results from other groups. The O'Neil & Bothun data in 
contrast to the Freeman result shows a flat distribution over 



9 CONCLUSIONS 

We have introduced the bivariate brightness distribution as 
a means by which the effect of surface brightness selection 
biases in large galaxy catalogues can be investigated. By cor- 
recting for the light below the isophote and including a sur- 
face brightness dependent Malmquist correction we find that 
the measurement of the luminosity density is increased by 
~ 37% over the traditional — method of evaluation. The 
majority (25%) of this increase comes from the isophotal 
correction with 10% due to incorporating a surface bright- 
ness dependent Malmquist correction and 1% due to the 
clustering correction. 

We have shown that our isophotal correction is suitable 
for all galaxy types and that isophotal magnitudes with- 
out correction severely underestimate the magnitudes for 
specific galaxy types. We also note that the redshift incom- 
pleteness suggests that predominantly low surface bright- 
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Figure 15. The solid line shows the final surface brightness func- 
tion, with all the surface brightness and clustering corrections. 
The points are offset towards the low surface brightness direction 
by 0.55 mag to crudely correct for the bulge. Also shown is the 
O'Neil & Bothun (2000) data, as crosses with error bars. The 
shaded region is where the visibility line shown in the previous 
diagrams starts to cross the data and so the data is no longer 
completely reliable. 



ness galaxies are being missed however we also show that 
these systems are predominantly low luminosity and hence 
contribute little to the overall luminosity density. This is in 
part due to the high completeness of the 2dFGRS (~ 91 %). 

We rule out the possibility of the 2dFGRS missing a 
significant population of luminous giant galaxies down to 
fl e = 24.5 mag arcsec -2 (or fi — 23.5 mag arcsec -2 ) and 
note that the contribution at surface brightness limits below 
(i e = 23 mag arcsec -2 is small and declining. Dwarf galaxies 
greatly outnumber the giants and the peak in the number 
density occurs at the low luminosity selection boundary. The 
implication is that the most numerous galaxy type lies at 
M > —16.0 mag. The galaxy population as a whole follows 
a luminosity-surface brightness relation (M oc (2.4±J;|)^t e ), 
similar but slightly steeper, to that seen in Virgo, in the 
Hubble Deep Field and in SdM galaxies. This relation pro- 
vides an additional constraint which galaxy formation mod- 
els must satisfy. 

We conclude that our measure of the galaxy contribu- 
tion to the luminosity density is robust and dominated by 
conventional giant galaxies with only a small (< 10%) con- 
tribution from dwarf (M > —17.5 mag) and/or low surface 
brightness giants (fi e > 23 mag arcsec -2 ) within the selec- 
tion boundaries(-24 < Mb < -15.5 mag, 18.0 < ^ e < 24.5 
mag arcsec -2 ). However we cannot rule out the possibility 
of a contribution from an independent population outside of 
our optical selection boundaries. 

Our measurement of the luminosity density is j = 2.49± 
0.20 x lO 8 /iiooL0Mpc -3 and using a typical cluster mass-to- 
light ratio leads to an estimate of the matter density of order 
Om ~ 0.24 in agreement with more robust measures. 

Finally we note that the bivariate brightness distri- 
bution offers a means of studying the galaxy population 
and luminosity-density as a function of environment and 
epoch fully inclusive of surface brightness selection biases. 
Future extensions to this work will include the measure- 
ment of BBDs and "population peaks" for individual spec- 



tral/morphological types, and as a function of redshift and 
environment. This step forward has only become possible 
because of the recent availability of large redshift survey 
databases such as that provided by the 2dFGRS. 
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APPENDIX A: TESTING THE ISOPHOTAL 
CORRECTIONS 

It is possible to model several different types of galaxy and 
compare the isophotal magnitude and the total magnitude 
as calculated in §4 with the "true" magnitude. The models 
are simple, assuming a face on circular galaxy, composed of 
a bulge with a de Vaucouleurs r 3 law (de Vaucouleurs 1948) 
and a disk with an exponential profile (see Eqn. 0). 



f^bulge 



/j, e + 1.40 + 8.325[(r/r- e ) 



1/4 



1] (Al) 



where r e is the half-light radius of the bulge, /x e is the ef- 
fective surface brightness of the bulge. Here we define it as 
the mean surface brightness within rj^. However, it is more 
useful to define galaxies in terms of their luminosities and 
bulge-to-disk ratios than their effective radii or disk scale- 
lengths. 

The magnitude of a galaxy and the bulge-disk ratio can 
be found in terms of the above parameters, by: 



M = -2.51og 10 [10-°' 4S + 10-°' 4D ] 
B 1 



B/T B + D 1 + !/(§) 
B — n e + 1.40 - 2.5 log 10 [7.227r r\] 
D = fi - 2.5 log 10 [n a 2 } 



(A2) 



where B is the magnitude of the bulge and D is the magni- 
tude of the disk. B/T is the bulge-to-total ratio. Given the 



t Note that the term "effective surface brightness" is sometimes 
defined as the surface brightness at r e , for Ellipticals the cor- 
rection between these definitions of effective surface brightness is 
1.40. 
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parameters M, B/T, fi e and fi , a galaxy's light profile is 
fully defined. 

To calculate the difference between the total and the 
isophotal magnitude it is necessary to find the fraction of 
light lost below the isophote. Since the intrinsic detection 
isophote varies with the redshift, this difference will be a 
function of redshift. For a variety of redshifts from z = 0.001 
to z = 0.201, the fraction of light under the isophote was 
calculated, by first converting the above magnitudes to ap- 
parent magnitudes, the intrinsic surface brightnesses to ap- 
parent surface brightnesses and then calculating the scale- 
lengths as above. The conversions from absolute to apparent 
properties are given in Eqn. [l7| and Eqn. [l8[ 

Using fiiim = 24.67mag arcsec -2 , the isophotal radii of 
the disk and bulge are calculated. 



/ f^lim M° \ 
Tiso d = ( ) Of 

' v 10g6 I 

_ f fMim - (Me + 1-40) n 4 
Tiso,b — { g ggej ' J- J ?"e 



(A3) 



The fraction of light above the isophote is then calcu- 
lated using the equation below. 



2/3-1 

/=i- v V e_9 

^ — ' n! 

n = 



(A4) 



where j3 is the de Vaucouleur's parameter, which is 1 for a 
disk and 4 for a bulge, g — 7.67( — ) 1,/4 in bulges and g = — 
in disks. The isophotal magnitude and isophotal radius of 
the galaxy can now be calculated. 



—21. and the values of Mj s0 and M corr are calculated at 
z — 0.12. Adjusting the value of Mtrue does not affect the 
changes in magnitude at any particular redshift, but does 
effect the maximum redshift that the galaxy can be seen at. 
The differences between Mi so & A/ true and M corr & M tr ue 
both increase with redshift. In all cases apart from the Ellip- 
tical galaxy, the value of M corr is a substantial improvement 
over Mi SO and remain below the photometric error for all 
types out to z = 0.12. 



APPENDIX B: VISIBILITY THEORY 

The equations below are reproduced from Phillipps, Davies 
& Disney (1990) and Disney & Phillipps (1983). These are 
used to calculate the volume over which a galaxy of abso- 
lute magnitude M, and central surface brightness \i can be 
seen. The theory determines the maximum distance that a 
galaxy can be seen to, using two constraints: the apparent 
magnitude that the galaxy would have, and the apparent 
size that the galaxy would have. The first constraint sets a 
limit on the luminosity distance to the galaxy which is the 
distance that a galaxy is at when it becomes too faint to be 
seen. 

dx = [f(»Um ~ Mo)]' IQ^^-M-K-K^ M P C 

(Bl) 

where f(fiu m — Mo) is the fraction of light above the 
isophotal detection threshold and is profile dependent. For 
a spiral disk with an exponential profile, the fraction of light 
is: 



m t30]b =mb- 2.51og 10 f b 
m iso ,d = m d — 2.51og 10 f d 
m lso = -2.51og 10 [10~ a4m " 

r iso ~ max ( r iso.b) *iso,d) 



+ 10 



(A5) 



Now that the isophotal magnitude, mi so, has been 
found it is possible to convert it back to an absolute mag- 
nitude, Miso- The isophotal magnitude and radius are fed 
now back through the equations in §4 and a value of m corr is 
calculated. This is converted to an absolute magnitude and 
Table 4 shows a comparison of M tr ue, Mi so and M corr at 
2 = 0.12 for the main Hubble Type galaxies and low surface 
brightness galaxies. The properties of the LSBGs (B/T, fi? 
and Ho) come from averaging the B-band data in Beijersber- 
gen, de Blok & van der Hulst (1999), the values of B/T for 
SO, Sa, Sb and Sc galaxies are taken from Kent (1985). The 
He values are taken from Fig. 5 of Kent (1985) for Ellipticals 
and Spirals, extrapolating where necessary (i.e. for Sb, take 
a value between that median for Sa-Sb and Sbc+) and sub- 
tracting 1.40 for the conversion from the surface brightness 
at r e to the effective surface brightness. For [i , we used the 
Freeman value for Spiral disks (Freeman, 1970). Irregular 
galaxies are usually placed beyond Spirals on the Hubble 
Sequence. They have either no bulge or a small bulge (van 
den Bergh, 1998; Smoker, Axon & Davies, 1998) and the 
majority are small Mb > —17. Ferguson & Binggeli (1994) 
show that all galaxies with Mb > — 16 can be fitted with 
exponential profiles. 

The absolute magnitude is kept constant at M trU e = 



fiHiim ~ Ho) = 1 - [1 + 0.4 ln(10)Gu« 



^ 10 [-0.4( mm -Mo)] 

(B2) 



The values of nu m and \i have to both be absolute or 
both be apparent, for the relation above to be true. As the 
known value of \iu m is the apparent value and the known 
value of fi is the absolute value, a redshift dependent factor 
must be included in all calculations. 

jUiim - Mo -> ihim - Mo - 10 log 10 (l +z)- K(z) 

(B3) 

Thus the maximum distance has a redshift dependence. 
The luminosity distance is a function of redshift: 



d L (z) = ^[(l + z)-(l + z) +0 - 5 ] (B4) 

tio 

The maximum distance can be found numerically by, 
for instance, a Newton- Raphson iteration as 



di(z)-dh{z) =0 



(B5) 



at the maximum distance. 

The second constraint, the size limit is found by a sim- 
ilar methodology. The formulation of the size limit is: 



d 2 = Cg(m im - Mo) 10 [(L2(fl °- M)1 /Bum Mpc 



(B6) 
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where C is a profile dependent constant. g(fiu m — 
is the isophotal radius in scale lengths. 9u m is the minimum 
apparent diameter. For a spiral disk with an exponential 
profile: 



C = 4= 10" 5 



2tt 



(B7) 



g^um - Mo) = 0.41n(10)Di«£ - Mo - 10 log 10 (l + 2) - K(z)} 

(B8) 

In this case di is an angular-diameter distance, not a 
luminosity distance. 



dAW = ^t(i+*r i -(i+*r 1 - 5 ] (B9) 



d 2 (z) - d A (z) = 



(BIO) 



Once the redshi fts Z \ and z%, which are the solutions 
of Eqn B5 and Eqn BIO , have been found, the maximum 



redshift which the galaxy can be seen to is the minimum of 
Zi and Z2- 



= min(zi, Z2) 



(Bll) 



For the 2dFGRS, the parameters above are: fiu m = 
24.67mag arcsec -2 , mu m = 19.45mag, and 9u m = 7.2". 
Zmax < 0.12, is another limit imposed on this survey. 

There are also limits on the minimum redshift. These 
are also caused by selection effects in the survey. One limit 
comes from the maximum apparent brightness of galaxies in 
the sample, which is m = 14.00. The minimum redshift, 23 
is calculated just as zi . There is a also a maximum radius 
for galaxies, which is 9 max = 200". This comes from the sky 
subtraction process. This leads to 24, calculated in the same 
way as 22- 



= max(z3, Z4) 



(B12) 



Zmin > 0.015. This limit is to prevent peculiar velocities 
dominating over the Hubble flow velocity. 



V(M,ii,) = 



a cd r 



H (l+z) 



rdz 



(B13) 



where a is the area on the sky in steradians, c is the speed of 
light and H is Hubble's constant. The visibility V(M,[i ) 
represents the volume over which a spiral disk galaxy of 
absolute magnitude M and central surface brightness /i D can 
be observed. The central surface brightness can be converted 
to effective surface brightness using the following formula. 



He = Mo + 1.124 



(B14) 
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Table CI. A comparison of the local luminosity density from recent magnitude-limited redshift surveys. 



* 2 Survey 


Reference 


M* 






a 


is/10 8 /iiooL Mpc" 3 


SSRS2 


Marzkc et al. (1998) 


-19.43 


1.28x10 


-2 


-1.12 


1.28 


Durham/UKST 


Ratcliffc et al. (1998) 


-19.68 


1.7x10" 


-2 


-1.04 


2.02 


ESP 


Zucca et al. (1997) 


-19.61 


2.0x10" 


-2 


-1.22 


2.58 


LCRS* 


Lin et al. (1996) 


-19.19 


1.9x10" 


-2 


-0.70 


1.26 


EEP 


Efstathiou et al. (1988) 


-19.68 


1.56x10 


-2 


-1.07 


1.89 


Stromlo/APM 


Lovcday et al. (1995) 


-19.50 


1.40x10 


-2 


-0.97 


1.35 


Autofib 


Ellis et al. (1996) 


-19.20 


2.6x10- 


-2 


-1.09 


2.05 


CfA* 2 


Marzkc et al. (1994) 


-19.15 


2.4x10" 


-2 


-1.00 


1.71 



* The LCRS used a Gunn r filter. The value of M* has been converted to Mb using < bj — R >o= 1.1 for the Johnson B-band (Lin et 
al. 1996). * 2 The CfA used Zwicky Magnitudes. The value of M* has been converted to Mb using < bj — Mz >o= —0.35 for the 
Johnson B-band and <fi* has been reduced by 60% (Gaztanaga & Dalton (2000). 



Table C2. A table of the number density as a function of M and fie for the 2dFGRS data. The values are in units of 1.0 X 10 _4 galaxics 
Mpc - 3 . The dashes are in bins where visibility theory says that the volume sampled, V{M,)i) < 10 4 Mpc 3 . The bold numbers are in 
bins where there are more than 25 galaxies. The italic numbers are in bins where there are less than 25 galaxies and visibility theory is 
used to calculate the number density. Where this becomes less than 1 galaxy in 10 4 Mpc 3 , the limit < 1 is used. 



Bin 


20.1 


20.6 


21.1 


21.6 


22.1 


22.6 


23.1 


23.6 


24.1 


24.6 


-24.0 










< 1 


< 1 


< 1 


< 1 


< 1 


< 


1 


< 1 


-23.5 


< 


1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 


1 


< 1 


-23.0 


< 


1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 


1 


< 1 


-22.5 


< 


1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 


1 


< 1 


-22.0 


< 


1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 1 


< 


1 


< 1 


-21.5 


< 


1 


< 1 


< 1 


1.4±0.4 


< 1 


< 1 


< 1 


< 1 


< 


1 


< 1 


-21.0 


< 


1 


< 1 


4.2±0.6 


6.1±0.8 


3.3±0.6 


< 1 


< 1 


< 1 


< 


1 


< 1 


-20.5 


< 


1 


1.9±0.4 


8.8±0.9 


18±1.3 


12±1.1 


3.0±0.5 


< 1 


< 1 


< 


1 


< 1 


-20.0 


< 


1 


2.9±0.7 


14±1.2 


42±2.0 


33±1.8 


9.0±0.9 


< 1 


< 1 


< 


1 


< 1 


-19.5 


< 


1 


< 1 


15±1.2 


59±2.4 


65±2.5 


26±1.6 


5.7±0.8 


1.3 


< 


1 


< 1 


-19.0 


< 


1 


1.0 


13±1.4 


58±2.4 


98±3.1 


52±2.4 


14±1.4 


4.2±1.0 


< 


1 


< 1 


-18.5 


< 


1 


< 1 


11±1.7 


61±3.5 


111±4.4 


77±3.7 


35±2.8 


11±1.9 


3. 





< 1 


-18.0 


< 


1 


< 1 


15±4.1 


43±3.5 


110±5.8 


113±6.0 


55±4.6 


21±3.4 


2. 


7 


< 1 


-17.5 


< 


1 


< 1 


2.8 


31±4.4 


96±7.4 


126±8.6 


62±6.7 


35±6.9 


7. 


1 


< 1 


-17.0 


2. 


5 


1.6 


2.2 


39±8.9 


113±13 


136±14 


103±14 


44±10 


20 


< 1 


-16.5 






< 1 


2.3 


27 


101±16 


151±19 


116±17 


66±14 


30 




-16.0 










19 


50±14 


171±27 


130±24 


72 








-15.5 
-15.0 














126±34 


102 










-14.5 



























Table C3. A summary of the luminosity densities and upper limits on Qm found from the 2dFGRS data. 



Parameters jr7l° 8 ^-^Q M P c 3 

HUm = 24.67 2.49 ±0.20 0.24 ± 0.05 

Hum = 24.97 2.33 ±0.17 0.23 ± 0.05 

H lim = 24.37 3.00 ±0.26 0.29 ± 0.06 

Simple LF 1.82 ± 0.09 0.18 ±0.04 
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Table C4. A summary of the errors in the isophotal and corrected magnitudes. 



Hubble Type 


Mtrue 


B/T 


Me 


Mo 


M iso 




E 


-21.00 


1.00 


20.5 




-21.00 


-21.07 


SO 


-21.00 


0.65 


19.2 


21.7 


-20.83 


-20.89 


Sa 


-21.00 


0.50 


19.6 


21.7 


-20.75 


-20.86 


Sb 


-21.00 


0.30 


20.4 


21.7 


-20.64 


-20.86 


Sc 


-21.00 


0.15 


21.0 


21.7 


-20.54 


-20.89 


Sd 


-21.00 


0.10 


21.2 


21.7 


-20.51 


-20.92 


LSBG 


-21.00 


0.12 


27.7 


23.0 


-19.08 


-20.84 


Irr 


-21.00 


0.00 




22.7 


-19.67 


-21.00 
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